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Abstract 

The transverse polarization of a quark is a degree of freedom that is not taken 
into account in the most commonly used Monte Carlo generators. For the case 
e~^e~ hadrons I show that it is possible to use these generators to simulate 
processes where the parent quark and antiquark are transversely polarized and the 
fragmentation process is affected by this polarization. The key point is that it is 
possible to obtain this without touching the generator code at all. One only works on 
the parton-level and hadron-level outputs that the Monte Carlo code has produced, 
modifying them in a correlated way. A group of techniques is presented to obtain 
this, matching the most obvious needs of a user (in particular: reproducing a pre- 
assigned final distribution). As an example these methods are applied to modify 
Pythia-generated events to obtain a nonzero Collins function and a consequent 
cos{^i + $2)— asymmetry of pion pairs. 
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1 Introduction 

A great effort covering several decades lias been devoted to developing full- 
purpose MonteCarlo codes (e.g. [T]|2|3] . see [1] for a more general review) for 
high-energy hard processes. However, the most known codes do not allow for 
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simulation of effects like azimuthal asymmetries in SIDIS, in e~^e~ — )■ hadrons, 
in Drell-Yan processes etc, i.e. those processes where the transverse spin of the 
quark has a relevance (for an overview of this field, see the proceedings [S], 
and the reviews [6] and [7]). 

Although it is possible to modify a Monte Carlo generator in such a way to 
produce such asymmetries, the perspective of touching the core of such com- 
plicate codes is surely a nightmare for most of the people that are potentially 
interested in. What I present here is a set of techniques aimed at introducing 
(transverse) spin physics in a generator without touching it at all. Of course 
one needs writing patches of code, but these use the parton level and hadron 
level outputs of the generatoiQ as an input, and are completely independent 
(small) programs. 

In this work, I will consider e~^e~ — ?• hadrons. I will start from events generated 
by Pythia-8 fL\. I need to modify Pythia outputs in such a way that: 

(i) In the hard vertex e~^e~ — ?• qq, transverse polarizations are added to the 
quark and to the antiquark, and the correlated distribution of spins and mo- 
menta is coherent with the known matrix elements. 

(ii) The momenta of the final hadrons in each hemisphere (quark hemisphere 
and antiquark hemisphere) are modified in a way that is correlated with the 
transverse spin of the parent quark. The modification must be under our full 
control, i.e. we should be able to obtain exactly what we want to obtain. This 
may mean a fragmentation function with pre-assigned form, or the implemen- 
tation of a model. 

I will describe some techniques to implement the previous requirements. As 
an example, I will apply them to modify pion final momenta so to have them 
distributed according with the sum of an unpolarized and a Collins fragmenta- 
tion function[8]. Although here I just want to show an example of application 
of these simulation techniques, the chosen case is of special interest, since the 
Collins function has a relevant role in the extraction of information on the 
transverse polarization of the nucleon[9], there are models for it [T0llll|12j . 
and it appears in asymmetries measured in SIDIS [T3|14|15j and e -f- e~ — >■ 
hadrons [T8] . 

As a consequence of a nonzero Collins function, the correlated distribution of 
pions detected in opposite hemispheres is expected to present a cos($i -|- $2) 
asymmetry in the sum of the azimuthal angles of the pions [TB] (see section 
VIII of [T7] for details). This will be confirmed by an analysis of the events 
generated by Pythia and modified as suggested here. 

^ The most known MC generators show, in their output, the momenta of all the 
particles produced in the intermediate stages of a complex event. 
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I will not try to reproduce the detailed physical outputs of the recent mea- 
surement of this quantity at Belle p^p9] . because this is just an exercise and 
the aim of this work is to present a more general group of techniques. Once 
the individual hadron tracks have been modified in a physically motivated 
way, other azimuthal asymmetries could be generated in a set of simulated 
events, like the cos(2</))— asymmetry (see [17] and [E]), or dihadron-dihadron 
correlations [20|2T] . or even asymmetries associated with a larger number of 
particles |22]. 

Let me name NPMC the "non polarized" MonteCarlo code whose outputs 
have to be modified by the external patches. 

The two main steps, corresponding to previous (i) and (ii) are: 

Step (1): Take a qq pair produced in the hard vertex by the NPMC and "stick" 
a pair of reciprocally independent random transverse spins, in such a way that 
the correlated spin-momentum distribution agrees with the polarized quark- 
lepton squared matrix element. Transverse spins are assumed to be classical 
fixed-length vectors with one degree of freedom (the angle in the plane that 
is normal to the quark momentum). 

This step is the critical point of the method. It needs to be demonstrated that 
it is feasible. One thing is to sort both momenta and spins according with a 
joint distribution, and another thing is first sorting momenta according with 
the spin-averaged distribution (that is done by the NPMC), and next dividing 
the sorted events into fairly distributed spin subsets (that is done by us). 
Section 2 is devoted to this. 

Step (2): Modify the (final or intermediate) hadronic momenta so to repro- 
duce the effect of an assigned quark-spin-dependent fragmentation function. 
Alternatively, one could like to implement a physical model that is behind this 
fragmentation function. 

Here the underlying assumption is that azimuthal effects are a small distor- 
tion of a final particle distribution that is mainly determined by the physics 
implemented in the NPMC. There are two classes of techniques that may 
be exploited: "distortion" , and "filtering" techniques. In distortion techniques 
the individual particle properties in a given event are modified. These tech- 
niques exploit each event produced by the NPMC. In filtering techniques only 
a subset of the events produced by the NPMC is accepted. 

Filtering may be expensive, when many NPMC-events are needed to obtain 
one final event. It may be necessary when we think that the additional physical 
processes affect quantities like the final pion multiplicity, or the general struc- 
ture of the event. If we do not expect this to be the case, distortion techniques 
are preferable. 
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In section 3, 1 show the effect of two possible distortion techniques in producing 
a nonzero Collins function for the pions. 



A cos($i + $2)— asymmetry involving pions of the opposite hemispheres is a 
synthesis of the presence of Collins functions on both sides, and of the spin 
correlations between the quark and antiquark produced in the hard vertex. 
So, if things in the steps (1) and (2) have been properly performed it must 
be present. Section 4 is devoted to show that this is the case, i.e. that a 
cos($i + $2)— asymmetry is present in a set of Pythia-events modified by the 
presented techniques. 



2 Step 1: Hard vertex and spin 



In most MontcCarlo codes for high-energy physics, the starting point is the 
hard scattering process at parton level. Once a hard scattering event has been 
sorted according with some probability, both later and previous cascading 
processes are generated. 



Method 1 




Fig. 1. Event distribution as a function of the quark- lepton angle 9 and of the 
components of spins of the quark and antiquark along the lepton axis. 

In the case of e+e~ — )■ hadrons, the (parton level) cross section for the hard 
scattering process e+e~ — )■ qq with transversely-polarized quarks, and with 
the leptons on the z— axis, is (see Appendix): 
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Fig. 2. Event distribution as a function of the quark-lepton angle and of the 
components of spins of the quark and antiquark on the plane normal to the lepton 
axis (see the next figure for the individual contributions of x and y terms). 




Fig. 3. Event distribution as a function of the quark-lepton angle and of the 
components of spins of the quark and antiquark along the x— axis (the corresponding 
y— distribution is equal). 



W^oi{N,sJ) = ^|[1 + cos\e)]{\ + s,-s,) 



- sin^{9){sxS^ + SySy) (1) 
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Fig. 4. Upper panel: s^s^— distribution integrated over the cos{6)— range (—0.5, 0.5). 
Lower panel: the cos (0)— distribution, summed over spins. The continuous line is 
obtained with Method 1, the dotted line with Method 2. 

where iV is a versor of the quark momentum (— is the corresponding one 
for the antiquark), 9 is the angle between the lepton and quark axes {cos{9) 
= N ■ z) and s, s are transverse polarization vectors. "Transverse" means 
"transverse to the (anti)quark 3-momenta". Since the quark momentum is 
not aligned with the electron momentum, the transverse spin has a nonzero 
Sz component, i.e. a component along the electron-positron beam axis. This 
must not be confused with the helicity or with the longitudinal spin of the 
quark. 

This cross section was already reported in [23J where it was applied to the 
Drell-Yan process in the context of the Panda experiment |24j . The two cross 
sections present a different form because they are expressed in different refer- 
ence frames, but they are the same. 

If the polarization components are averaged away, we get the known result 



WnopoiiN) oc 1 + cos\e). (2) 

The most known and used MonteCarlo codes implement this equation for 
the lepton-quark vertex, adding further 0(a") hard processes in the partonic 
showers. 
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Fig. 5. Upper panel: SxSx + SySj^— distribution integrated over the cos(0)— range 
(—0.5,0.5). Lower panel: the cos(0— )distribution, for —0.2 < SxSx + SySy < 0.2. 
The continuous line is obtained with Method 1, the dotted line with Method 2. 



If we had to write from the very beginning a complete MonteCarlo generator 
including transverse polarizations of the quark and antiquark, we could imag- 
ine two ways of sorting momenta and spins of the quark and of the antiquark 
(generically: "quarks" ) . 

Method 1) The momenta and the spins of the quarks are jointly sorted ac- 
cording with the probability law eqU] 

Method 2) First, the momenta of the quarks are sorted accordingly with eqJH 
and next the polarizations are sorted in such a way to get the distribution [T] 
for the joint set of variables N,s,s. 

The former one is the right one in general (if correctly implemented). The 
latter method allows for splitting the generation process into two steps, one 
of which may be performed by an NPMC and the other one by our additional 
code patches. I will show that in the case interesting us it produces the same 
results of method 1. 

To apply Method 2 here I take a quark-antiquark pair whose momenta have 
been sorted with probability Wnopoi{N), and sort their spins according with 
the distribution 
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WnopoliN) 

subject to the conditions 



s-N = s-N = 0, 
\s\ = 1, 1^ = 1. (4) 

In eqjni in general, D requires the denominator Wnopoi- In this pecuhar case it 
may be removed since we sort a variable that is absent in this denominator. 

For method 1 the implementation is more straightforward. Events are sorted 
according with the probability distribution eq. [H subject to the condition HI 

In the set of figures [TH5] the distributions of some variables and of their corre- 
lations are reported, calculated with both Method 1 and 2. These show that 
Method 1 and 2 produce the same results. 

Since the full distribution shows a strong correlation between 6 and SzSz on one 
side, and between 9 and s^Sx + SySy on the other side, the scatter plot of these 
pairs of variables is shown in figs. [T]|2J These figures report the distribution 
of 100,000 events in the (cos6', A)— plane, where A is s^s^ (figH]), s^Sx + SySy 
(figJlD, SxSx (figH. 

The other figures show slices/integrals of the previous distributions. In figJH 
(upper panel) I report the s^s^-distribution integrated over — 0.5 < cos{9) < 
0.5. This region is relevant for the Collins effect that is suppressed for 1005(6*) | 
near 1. In the lower panel the spin-integrated cos(^)— distribution is reported. 
It reproduces the 1 + €03"^ (9) shape within 5 %. In figE] (upper panel) I report 
the distribution of SxSx + SySy integrated over —0.5 < cos9) < 0.5. In the 
lower panel the distribution is the one of cos{9) for the integration range —0.2 
< SxSx + SySy < 0.2. In the last case, the larger (fluctuating) discrepancies 
between the distributions obtained by the two methods have statistical origin, 
and are due to the relatively small number of events in each bin. 

Concluding this part, I may claim that Methods 1 and 2 give similar results, 
so it is licit to apply the polarization stage of method 2 to events sorted by 
some independent generator that does not include spins. 



3 Step 2: Reasonable choices for the distortion methods 



Now we need to modify the transverse momenta of the final hadrons in such 
a way to reproduce an assigned distribution, that is correlated with the quark 
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polarization. In an attempt to be as comprehensive as possible, I propose and 
test two classes of methods, that I name "product" and "convolution" . These 
should fit the most obvious requirements. 

In this section, I apply these methods to a set of 2-momenta that are gaussian- 
distributed in the xy— plane with center of the distribution in the origin. They 
are supposed to be the transverse momenta of a set of final hadrons, originating 
from a parent quark that is directed along +z and has polarization +x. After 
testing the distorting methods on this simplified set, in the next section they 
will be applied to a set of Pythia events. 

Let Fq{Kt) be the undistorted momentum distribution, that depends on Kt 
via \Kt\ only. It depends on the longitudinal fraction Z (not explicitly re- 
ported) and is proportional to the fragmentation function Di{q — )■ h; Z, Kt), 
where q is one or a group of quark flavors and h is one or a group of hadron 
species. Fq is normalized to 1, while Di is normalized to the total hadron 
multiplicity in the subset of events q ^ h. 

Let G{Kt) be a distorting factor, that introduces azimuthal asymmetries in 
the xy— plane. It may depend on Z, but I do not write this explicitly. 

Product techniques: the final distribution has the form 

F{Kt) = Fo{Kt) G{Kt). (5) 

Convolution techniques: the final distribution has the form 

F{Kt) = J dkTdqrSiKT - kr - qx) Fq{Kt) G^qr). (6) 
or, more in general. 



F{Kt) = j dltTdqT^^Kr — kx — qr) 

Fo{KT)G{KT,qT). (7) 

The latter form is not strictly a convolution, but I will use this name anyway. 
Sometimes, a starting model suggests a form like in eqJHl but some constraint 
on KT + qr (when sorting (Jt) may remove the full independence of G{qT) from 
Kt. 

The product form allows an easy implementation of parametrizations of the 
form 

H{z,Kt) 



D,{z.Kt) 1 + 



Di{z,Kt) 



9 



The convolution form is more flexible to implement physical models, since it 
treats the flnal hadron momentum as a sum of momenta with independent 
physical origin. 



3. 1 Implementation: the product case 

To implement the product form I have chosen an algorithm belonging to the 
Metropolis-Hastings family [25|26] : 

1) I start with an "undistorted" event Kt- 

2) A shift qT is sorted (flat distribution) inside a circle of radius Rq. 

3) The "shift probability" is 



If this quantity is > 1, the step to the new point + Qt is performed. If it 
is not, an accept/reject procedure is set. As a consequence, with probability 
P{Kt, Kt + Qt) the step is performed, and with probability 1 — P{Kt, KT + qr) 
the step is not performed. 

4) As a result, we have a new point that either coincides with Kt or with 
Kt + Qt- 

5) The sequence 2-4 is repeated starting with the value K^ instead of Kt- A 
new point K^ is selected. 

6) The sequence 2-4 may be performed times. 
Comments: 

For A^ = 1 and for a small admitted displacement qr the flnal distribution is 
similar to the starting one Fq{Kt). 

For A^ = 10 the flnal distribution coincides with F{Kt) whichever the starting 
distribution was (e.g. one may choose Kt = flxed). So the joint choice of A^ 
and Fq{Kt) must be clever, in those cases where one does not know precisely 
the distribution of the undistorted events. 
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3.2 Implementation: the convolution case 



The implementation of the convolution method is more straightforward. In 
general terms, to get a distribution of the form h[X) = J dx J dyf{x)g{y)6{X — 
X — y) the steps are 

1) Sort X according with /(x). 

2) Sort y according with giy). 

3) Sum X and y. 

This may work to produce both a distribution of the form eq. [6l and one of 
the form eq. [71 



3.3 The examples of figlB 



As an example I have chosen, both for the undistorted distribution, and for the 
distortion factors, some shapes that simplify much the computational work (in 
particular, imposing that the distorting factors are zero in some parts of the 
phase space). Apart for this, they do not present any special lack of generality. 

In the examples of figEl as an undistorted and axially symmetric distribution 
Fo{Kt) I have chosen the gaussian 

Fo{Kt) = exp[-{Kl + Kl)/{2kl)] (10) 

with A;i = 0.5 GeV/c. 

For the product case I have applied eq. [5] with 

G{Kt) = 
1 + aKyexp[-{Kl + K'y)/{2kl)]^ 

when this expression is > 0, and 

(11) 

G{Kt) = 0. 

(12) 

when it is negative. 
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For the example in figj6]a = 0.75 and A;2 = 4 GeV/c >> ki. At any step, the 
shift qx is sorted with Qx = and \qy\ < 1 GeV/c. I have used = 10 steps, 
but aheady with 5 steps the final results converges towards the final required 
distribution. 



For the convolution case I have used 



a Qy exp[-{ql + g^)/(2A;^)] for Qy 



> 



Giqr) = for qy < 



0, 
0. 



(13) 



with a = 0.5 and = 0.1 GeV/c. These parameters are chosen so that the 
convolution produces a shift that is similar to the product case. The choices ^2 
» ki (product case), and k^ «ki (convolution case) lead to a similar result, 
i.e. avoiding that the final distribution is much broader than the starting one. 
This permits us to appreciate the shifting effect, that in fig. E] is highlighted 
by vertical lines showing the average point of each distribution. 



K„ distributions with/without Collins effect 
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Fig. 6. Distribution of 40,000 events obtained by sorting a gaussian distribution 
w.r.t. k-T (fine-dotted curve), and next by shifting events according to either the 
convolution method (continuous) or the product method (coarse-dotted). Vertical 
lines show the average point of each distribution (the two averages of the shifted 
events practically overlap). 
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4 Example: A cos($i + $2)— asymmetry in the complete event 



As an example, the two described track-distortion methods, with the same 
parameters as in figE] are applied to produce a cos(<l>i + $2) nonzero term 
in the distribution of pion pairs generated by Pythia in the conditions of the 
Belle experiment in the sub-66 energy range, at Q = 10.0 GeV (Belle has 
performed measurements in this range, and at a slightly higher Q at the bb 
threshold [19]). 

To apply the previous methods to events produced by Pythia some further 
complications are needed. I do not give details on these points since they just 
require standard operations like frame rotations and vector projections. What 
is done is (1) the distortion of the hadron transverse momentum is produced 
in a frame where the parent quark (or antiquark) momentum and spin are 
along the z and (this reproduces the situation analyzed in the previous 

section), (2) for the specific data analysis (i.e. extraction of the cos($i + 
$2)— distribution) each resulting hadron momentum is transferred to a frame 
where the lepton-quark scattering plane coincides with the a coordinate plane. 

The Pythia-generated events are selected by the additional condition Thrust 
> 0.8. This excludes three-jet events, and in the case of Q slightly over 10 
GeV (i.e. over the value used here) it would exclude events of 66— kind (see 
the discussion in [H]). This cutoff is important, because the momentum-spin 
correlation calculated in the Appendix of the present work refers to "light" 
quarks, i.e. fermions for which it is licit to assume helicity conservation in the 
vector-fermion vertex. At Q = 10 GeV we have to include events starting from 
uu, dd ss and cc pairs. In the collision cm. frame, the quark energy is 5 GeV. 
The charm mass is ~ 1.27 GeV/c^. For this quark the first order correction 
to the UR relation E = P is m'^/2E ^ 160 MeV << E, so the helicity 
nonconserving terms have small relevance. If gluon radiation processes with a 
large Pt were included in the hard photon-quark vertex, the spin-momentum 
correlations would not be exactly as in eq. [T] (e.g. the quark momentum would 
be modified by a non-collinear gluon radiation, and the large quark virtuality 
would affect helicity conservation in the photon vertex). All these processes 
are excluded by the request Thrust > 0.8. 

According with e.g. [H] in presence of nonzero Collins functions on both sides 
one expects a cross section of the form 



The angle 9 is the lepton-quark polar angle. $1,2 are the angles of the pion 



da 



oc 1 + cos^{9) D{zi)D{z2) 



dzidz2dVtidVt2 



•2 



+ sm(^)cos($i + ^2)C{z^)C{z2). 



(14) 
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transverse momenta w.r.t. the lepton-quark scattering plane. D{zi^2) and C{zi^2) 
are the unpolarized and Colhns fragmentation functions with opposite hemi- 
sphere pions with longitudinal fractions zi and Z2- 

Since this is just an example of application, I have not attempted reproducing 
the recently measured values for this quantity at Belle |18p9] . Rather, I have 
reported the output corresponding to the distribution shifts reported in figlHl 
including the "no-shift" case to have an estimator of the fake asymmetries due 
to the statistics and to the cutoffs. 

((|)^+(|)p-event distribution (convoiution methiod) 

66001 




Fig. 7. Distribution of the angle $i -|- $2; where the angles $j are the az- 
imuthal angles (see text) of two opposite-hemisphere pions. The individual pion 
tracks are deflected according with the convolution method. The fitting curve is oc 
1 + ^ cos($i + ^>2), with A = 0.04 ± 0.004. 

The three figures [71 El [9] show the distribution of 135,000 pion pairs vs cos{^i + 
$2) in Belle conditions. Each event has Thrust > 0.8, and each individual pion 
has z > 0.1, and cos{a > 0.8, where a is the angle w.r.t. the quark axis (in an 
attempt of precisely fitting the experimental data, one should include further 
cuts and use the Thrust axis instead of the quark axis). A "pair" is composed 
of two pions (regardless of their charge) belonging to different hemispheres. 

The produced histograms have been fitted by the Root-Migrad package [27] 
with curves of the form B[l + A cos($i -|- $2)]- The data are divided into 40 
equal-range bins, and B = 6130 ± 20 is the average number of events per 
bin. A has size 3-4 %, compared to a statistical error 0.4 % and to a zero- 
asymmetry value 0.6 % (extracted from the last figure, where no azimuthal 
effect is included). 

The chosen form of the spin- dependent fragmentation functions is aimed at 
simplifying the implementation of the presented examples. I have not con- 
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Fig. 8. Distribution of the angle $1 + $2) as in figiT) The individual pion tracks 
are here deflected according with the product method. The fitting curve is oc 
1 + A cos($i + $2), with A = 0.03 ± 0.004. 



((|)^+(|)^)-event distribution (no shift) 
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Fig. 9. Distribution of the angle $1 + ^2, as in figiTl The individual pion tracks are 
not deflected at all, i.e. they are the original trajectories as simulated by Pythia. 
The fitting curve is cx 1 + A cos{<^i + ^>2), with A = -0.006 ± 0.004. 



sidered a dependence of the azimuthal effects on the longitudinal fractions. I 
have not considered the flavor dependence of the spin-dependent fragmenta- 
tion functions. I have only considered spin effects on pions. To include different 
functional forms, a dependence on z, flavor differentiation, azimuthal effects 
on other hadron species (like kaons), would be a matter of more program lines, 
but would not issue any further challenge. 
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The same procedure could have been apphed to an intermediate state hadron 
that decayed into two hadrons hke e.g. a p. Then the two decay hadrons would 
present an extra fragmentation function of the kind (see ^28j, t29j). On 
the other side, to produce a function of the kind H'^ ([2S], [2S], in e"'"e^ — )■ 
hadrons it has been measured by [30]) modifications should be applied to 
the relative momentum of the two final hadrons. These are straightforward 
generalizations of the work presented here. Another possible application is the 
production of polarized A/A pi]. In this case one would sort a spin for this 
hadron correlated with the quark spin, and sort new momenta for the decay 
products according with the assigned spin. 



5 Conclusions 

I have started this work from events of the class e"'"e~ — )■ hadrons, generated 
by an ordinary Monte Carlo code where the transverse polarizations of the 
quarks are not taken into account. 

A set of techniques have been presented to build code patches that do not 
touch the Monte Carlo generator itself, but rather modify its outputs, both 
at parton and hadron level. 

The first aim of these modifications is attributing a transverse spin to the 
quark and to the antiquark pair produced in the hard vertex, in such a way 
to get a physically sound distribution for the correlated set of quark momenta 
and spins produced in the hard vertex. To reach this, e"'"e~ events with pro- 
duction of polarized quarks have been generated according to a known and 
general method. Next, an alternative two-step method for generating similar 
events has been tested, where in the first step an unpolarized quark and an 
unpolarized antiquark are produced in the lepton annihilation, and in the sec- 
ond step both are polarized without touching the previously sorted momenta. I 
have shown that the two methods produce the same distributions. This means 
that it is licit to take an unpolarized quark- antiquark pair produced from an 
external generator and attribute a pair of polarizations to it according to the 
second step of the two-step method. 

The second aim is to slightly modify the transverse momenta of the produced 
hadrons, in a way that is related to the spin of the parent quark or antiquark, 
and controllable. "Controllable" may mean two alternative possibilities: Ei- 
ther that a given physical model is implemented, or that we know the form of 
the momentum distribution, or of the fragmentation function, that we want 
to obtain. Two classes of techniques have been presented, namely "product" 
and "convolution" techniques. The former group is suitable for producing final 
distributions according with fragmentation functions of the form D{1 + H/ D), 
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where D is the unpolarized fragmentation function and H a spin-dependent 
term, hke e.g. a Colhns function. The "convolution" techniques are more ap- 
propriate for those cases where a model predicts that the final moment is a 
convolution of two contributions, one due to the unpolarized physics, and the 
other one to polarization-related effects. 

As a test case, this has been applied to modification of a set of Pythia events, 
so to produce an azimuthal asymmetry that derives from the combined effect 
of the Collins functions of two opposite- produced pions. This was just a test 
case, the range of possible applications is quite large. 

As a final and due observation, I remark that the events produced by the 
methods described here, working on Pythia outputs, are not Pythia events. 
They are modifications of Pythia events. 
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6 Appendix: Polarized quark — unpolarized lepton contraction H'^'^{sy, Sy)L^ 

For the procedure described in section 2 we need the contraction 

q^W = H'^'L^, (15) 

of the quark-level hadronic and lepton tensors in the processes — )■ qq, 
where the final pair is constituted by a transverse-polarized quark and anti- 
quark. This will be calculated for purely massless quarks and leptons, in the 
center of mass frame of the reaction. 

EqJTSl gives, for assigned momenta of the electrons, the joint probability for 
sorting the momenta and the spins of the quark and the antiquark. In a lepton 
annihilation process, q has no probabilistic role since it is fixed. 

In ref. [23] (see the Appendix of that work) the q'^W factor was already calcu- 
lated for the reversed process qq — )■ /"*"/", aimed at the Drell-Yan application. 
The invariant result is the same in both cases, however relevant differences 
appear when one rewrites it in the reference frame where the transverse spins 
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need to be generated. The present case is simpler than the Drell-Yan one 
and requires no approximations, since the leptons are always on the 2;— axis, 
and the polarized quark and antiquark are exactly back-to-back. Because of 
this, the transverse spins of the quark and of the antiquark are orthogonal 
to the same axis. The 4-vector associated to the quark spin is orthogonal to 
the 4-momenta of the quark and of the antiquark. The former orthogonality 
is always true, the latter only for back-to-back pairs. 

I use the shortened notation for traces 

T[..] = \Tr[..]. (16) 

I use the definitions of the Berestevskij-Lifsits-Pitaevskij book [32] (better 
known as the 4th book of the Landau-Lifsitz Course in Theoretical Physics). 
As the only exception to this, I indicate k^'j^ with the widespread notation ^ 
instead of using k as was done in that book. 



6. 1 The case of unpolarized quarks 

If the quarks are unpolarized, we simply have 



t 
2 



- (17) 



and 

L^^ ^ {V^p"} - (18) 

where braces indicate symmetric dyadic product. The contraction of the two 
is faster using 

{kMiv^.f} = 2[{kp){kp) + mmi (19) 

[{k,,K} - {kk)g,,^ g^^ = 0. (20) 

We get 

2 [{kp){kp) + ikp){kp)]. (21) 
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I extract a factor q from each vector: 

P = {q/2)N^\ = (g/2)n^,... (22) 

In a center of mass frame of the partonic process, the 3-vectors A^, n etc are 
unitary vectors. I get 

Wunpoi = \ [{Nn){Nn) + {Nn){Nn)]. (23) 
o 

A more famihar way to write this may be 

Wunpoi =\[^ + cos\e)]. (24) 

where 9 is the angle between the lepton and quark directions in the partonic 
center of mass frame. 

6.2 polarized quarks 

In the case of polarized (anti) quark, we need to substitute 

^ ^ ^(1-7^^) (25) 
where is the polarization 4-vector for the quark, respecting the exact 4- 
dimensional constraint {sk) = 0. If each spin is exactly transverse to the 
corresponding momentum, then also s-k = 0. 

The relation between and the polarization in a rest frame a is, for massive 
particles. 

So = crL\k\/m, Sl = (JLE/m, St = (?r (26) 
and evidently it creates problems for m ^0, unless the longitudinal component 
is strictly zero. However, when the previous expressions are used to write the 
density matrix eql25] in terms of the rest frame polarizations, E/m— terms 
cancel and we may check that the density matrix is free from mass singularities: 

^(1-7^^) = ^[l-7'(±t^L + c?T-7T)]. (27) 
(± differentiates particles and antiparticles). 

Although it is not the aim of this particular work, it is useful to consider what 
would be the outcome for longitudinally polarized quarks, with helicities h 
and h: 



Wunpol Wf^fi = 

{1 - hh)Wunpoi oc {l-hh)[l + cos^e)]. (28) 
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This is a predictable and known result for massless fermions: pairs with the 
same helicity are suppressed since they have opposite longitudinal spins in the 
cm., i.e. total longitudinal spin zero. This cannot be transferred immediately 
to the transverse spin, since any helicity component is composed by 50 % 
± components along any chosen transverse axis. In addition, the 1 + cos'^{6) 
factor is a composition of Yoo and Yiq eigenfunctions of the orbital angular 
momentum along the 2;— axis, not along a transverse axis. So, while eqj28] 
could be guessed from the very beginning, a statement like "the transverse 
spins of the pair are mostly parallel" has not a solid a priori justification. 

The full trace for polarized quarks is 

^ r[^7^(l-75^)^7'^(l-7V)] (29) 

that excluding O ^(7^)^^ terms (these would lead to an antisymmetric tensor, 

that is useless when contracted with the symmetric tensor of the unpolarized 
leptons) reduces to 

= T[lj^-f^'l^-f''] - Tll^-f" P]- (31) 
Next I use the reduction formula 

T{abcdef) = g''^T{cdef)-g^^T{hdef) + 
+ g'"^T{hcef) - g^^Tihcdf) + g^^T{bcde) (32) 

(with the obvious notation T{ab...) = T[7"7''...]). For any of the remaining 
traces we have the better known relation 

T{abcd) = g^'^g^ - g^^g^d ^ g^dgbc ^33^ 

The composition of eqs. [32] and [33] produces 15 terms. Applied to the second 
term of eq. [311 two of these terms contain the products (sk) = and (sk) = 
and are dropped. 

In the specific case of the transverse spins for a back-to-back qq pair we also 
have {sk) = and {sk) = 0. So most terms are dropped and we are left with 



- (ss)] ■ HZ,oi - |K'^^}j (34) 

After contracting the hadron tensor with the (unpolarized) lepton tensor, the 
final result is 
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W,^s-r = \\\^ + cos\9)]{l + sj,) 

- sin^{9){s^s^ + SySy)j (35) 

Here I have used eqlMJ The appearance of nonzero s^, Sz terms should not 
confuse about the transverse nature of the considered spins. These terms ap- 
pear because the 2;— axis is parallel to the lepton momentum, not to the quark 
momentum. 

The SjSj— terms may have two different origins in the previous invariant equa- 
tions: from terms like (ss), and from terms like {sp){sp). When passing to 
3-dimensional components, both terms change sign. The former because {ss) 
— i- —{s-f). The latter term presents two such changes of sign, and a third due 
to the opposite space parts of and p^. 
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